Shape transition and propulsive force of an elastic rod rotating in a viscous fluid 



O 
O 

<N 

O 

Q 

m 



o 

■i— > 



O 

o 



> 

m 
o 

r-- 




Bian Qian, Thomas R. Powers, and Kenneth S. BreueiB 

Division of Engineering, Box D, Brown University, Providence, RI 02912, USA 

(Dated: February 2, 2008) 

The deformation of thin rods in a viscous liquid is central to the mechanics of motility in cells 
ranging from Escherichia coli to sperm. Here we use experiments and theory to study the shape 
transition of a flexible rod rotating in a viscous fluid driven either by constant torque or at constant 
speed. The rod is tilted relative to the rotation axis. At low applied torque, the rod bends gently 
and generates small propulsive force. At a critical torque, the rotation speed increases abruptly and 
the rod forms a helical shape with much greater propulsive force. We find good agreement between 
theory and experiment. 

PACS numbers: 87.16.Qp, 47.15.G-, 46.70.Hg 



Understanding how flagella and cilia work is a central 
aim of the field of cell motility. The problem may be split 
into two parts, both of which involve physics: the means 
of actuation, and the fluid-structure interaction. In this 
letter we consider the fluid-structure interaction for thin 
filaments in a viscous fluid. At micron scales, viscous ef- 
fects dominate inertia, and the fluid-structure interaction 
problem simplifies because the Stokes equations govern- 
ing the fluid motion are linear. Gray and Hancock used 
this linearity to develop a simple theory that successfully 
predicted the swimming speed of a sperm cell with a load- 
independent pattern of bending waves propagating along 
the flagellum [1]. Soon after, Machin considered the 
fluid-structure interaction [2[. He argued that the mo- 
tors must be distributed all along the length of the flag- 
ellum, since, for small amplitudes, a passive flexible rod 
waved at one end has an exponentially decaying envelop 
of deflection, whereas the amplitude of deflection in real 
flagellar bending waves increases slightly with distance 
from the head [2] . The shapes and propulsive forces of a 
passive rod actuated at one end have recently been exam- 
ined theoretically [3, [4J and experimentally [4| . Although 
sperm flagella are not passive, the results of @, 0, 0] are 
important for modeling real flagella since the modes that 
Machin found also enter models in which the flagellum is 
actuated along its entire length [5[. 

Rotating flagella are also common. For example, nodal 
cilia Q have an internal structure similar to that of sperm 
flagella. However, instead of beating in a plane like most 
sperm flagella, nodal cilia rotate along the surface of an 
imaginary cone. The flow set up by these flagella has 
been implicated in the formation of left-right asymmetry 
in developing embryos (see @ and references therein). 
Bacterial flagella provide another example. These flag- 
ella are helical, much thinner than eukaryotic flagella, 
and driven by a rotary motor embedded in the cell wall. 
Fluid-structure interactions are important for polymor- 
phic transformations in swimming bacteria [7] and the 
bundling of multiple flagella [§]. 

Complementary to the problem of understanding how 
biological flagella work is the problem of building an arti- 



FIG. 1: Orthogonal images of steady-state shapes of rotating 
rod with torque just below (a) and just above (b) the critical 
torque. The motor (not shown) is at the top, with rotation 
axis along z. Gravity points down. In (a) and (b), the left 
panel is the side view, and the right panel is the front view. 
The rod is marked with white dots for contrast. The axes in 
(b) are the same as in (a). The curved arrow in (b) denotes 
the sense of rotation of the rod. 



ficial microscopic flagellum-propelled swimmer. This feat 
was recently accomplished by Dreyfus et al., who used a 
rotating external magnetic field to generate propagating 
planar bending waves in a filament composed of a string 
of colloidal magnetic particles [9]. A major challenge in 
building an artificial microscopic swimmer is the means 
of actuation. Manghi et al. proposed a simple mech- 
anism in which a microscopic flexible rod rotates along 
the surface of an imaginary cone [l(J EH • Using numer- 
ical methods, they predicted that at a critical driving 
torque the rod will undergo a discontinuous transition to 
a helical shape with significant propulsive force, indepen- 
dent of the sense of rotation. In this letter, we present a 
macroscopic experimental realization of this system. We 
also present new theoretical results that complement the 
hydrodynamic calculations presented in ToL 111] . 

In our experiment, a servo motor rotates a flexible rod 
in highly viscous silicone oil. The rod is connected to 
the motor shaft with a kink such that the base of the 
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FIG. 2: Lumped parameter model consisting of two rigid 
links connected by a torsional spring (open circle). The top 
link is clamped. All drag is concentrated at the two filled 
circles. 



rod makes a fixed angle with the rotation axis (Fig. QJi). 
The motor may be operated either at constant speed, or 
at constant torque using a torque sensor with feedback. 
The range of torques explored was 0.5 to 8 mN-m, and the 
maximum rotation frequency was less than 0.3 Hz. The 
rod is a steel extension spring wrapped in Teflon™ tape; 
the tape stiffens the rod to minimize sagging at low speed. 
The diameter of the rod is a =2.5 mm, and the bending 
modulus is A w 3 x 10 _3 N/m 2 . Various rod lengths L 
from 210 mm to 290 mm were tested. The silicone oil 
has viscosity rj « 110N-m 2 /s and is held in a tank 420 
mm on each side. With these parameters, the Reynolds 
number Re = pvL/r] « (10 -2 , where p « 10 3 kg/m 3 is 
the fluid density, v ~ 10 _1 m/s is the typical velocity of 
the free end of the rod, and L « 10 _1 m. Front and side 
images of the steady-state three-dimensional shape of the 
rotating rod at each torque was captured using a single 
camera and a single mirror. The imaging system was 
carefully calibrated to account for perspective, achieving 
an accuracy of ±2 mm. 

At low torque, the rod bends slightly and the rotation 
speed is relatively slow (Fig. [T^l) . Just above a criti- 
cal torque, the rod adopts a helical shape and rotates 
much faster (Fig. [TJd). We first give a simple analysis of 
this shape transition using the lumped parameter model 
shown in Fig. [2j It is convenient to imagine prescrib- 
ing motor speed uo rather than motor torque. Since the 
Reynolds number is small, we take Re=0. Thus, we may 
work in the rod's rotating rest frame without introduc- 
ing fictitious forces. The flow in this frame at point r 
is uoz x r. The rod is modeled by two rigid links of 
unit length connected by a torsional spring. The link 
OP is fixed and represents the condition of fixed an- 
gle between the rotation axis z and the base of the 
rod in our experiment. The torsional spring represents 
bending resistance and is only sensitive to changes in 
the angle between the vectors OP and PQ. Assuming 
<C 1, the moment about P on PQ from the spring is 
M b « K(OP x PQ) « K(y, 20 - x, 0) , where K is the tor- 
sional spring constant, and (x, y, 2) w tq is the position 
of the point Q to leading order in 0. To find the steady- 
state position of Q, equate the moment on PQ due to the 
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FIG. 3: (Color online) Dimensionless motor torque MmL/A 
was measured as a function of dimensionless speed \ f° r L — 
250 mm (O) and L 290 mm (A) with angle = 26°. For 
L = 250 mm, speed was measured as a function of increasing 
torque (■) and decreasing torque (♦). Note the hysteresis. 

The linear ( ) and nonlinear ( — ) predictions are shown. 

The insets show examples of the steady state filament shapes 
in the low (left) and high (right) speed regimes. 

torsional spring to the moment on PQ due to the flow. 
Assuming all drag on PQ is concentrated at Q (Fig. [2]), 
the viscous moment about P is M v w — (uo(x, y, 0), where 
C is a resistance coefficient. Solving moment balance for x 
and y yields x = 20 /[I + ((uo/K) 2 } and y = x(uo/K. As uo 
increases from zero, the link PQ deflects and y increases, 
which causes Q to experience a viscous force in the nega- 
tive x direction. These forces push Q toward the rotation 
axis, and tend to cause the rod in our experiment to wrap 
around the z axis. As uo increases further, Q moves close 
to the rotation axis, and y begins to decrease. There is 
also some drag on the link OP, concentrated a distance 
d from O. The moment about O due to flow is 

Mo = ^ 2 ( rf2 + i + c^ 2 )- (1) 

For d 2 < 1/2, we find that the moment first increases 
with cj, then decreases as the link folds in toward the 
rotation axis where the flow is slow, and then increases 
again as the drag from the base link OP dominates. If 
Mo is plotted vs. oj, then we find an S-shaped curve, 
just as in our experiment (Fig. [3|), with discontinuous 
transitions in shape and speed as moment varies. 

We now turn to a quantitative analysis of the experi- 
ment. We will continue to prescribe uo rather than mo- 
tor torque M m in our calculation, and limit the analysis 
to steady-state shapes. Unlike Manghi et al. [ToL fllj. 
we disregard hydrodynamic interactions between distant 
parts of the rod and use resistive force theory to model 
the force per unit length f acting on the rod p], [l2[ : 

f = Cl(v - r s r s • v) + C||r s r s • v, (2) 
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where (± = 4^7?/ [log (L /a) + 1/2] and C|| = 
27rr]/[log(L/a) — 1/2] « C-l/2, r(s,£) is the position of 
the point on the rod centerline with arclength s at time 
t, v(s,t) is the velocity of the undisturbed flow relative 
to the velocity of the rod at s, and r s = dr/ds is the 
tangent vector to the rod centerline. There is also a hy- 
drodynamic torque per unit len gth distributed along the 
rod that tends to twist it [13, LLJ, LL5[ . However, the effects 
of this torque are smaller by a factor of (a/L) 2 relative to 
effects due to translation of the rod [13[ and will hence- 
forth be disregarded. The constitutive relation for the 
elastic rod is 

M = Ar s x r ss , (3) 

where M is the moment due to internal stresses exerted 
on the cross-section of the rod at s, and A is the bending 
modulus [16]. The shape of the rod is determined by 
force and moment balance, 

F s + f + f g = (4) 
M s + r s x F = 0, (5) 

where F(s) is the force due to internal stresses acting on 
the rod cross section at s, and f g = (/i ro d — ™ 2 p ii)#z 
is the buoyancy force per unit length due to the den- 
sity difference between the rod (linear density /i ro d = 
0.0478 kg/m with oil inside) and silicone oil (p n = 
970kg/m 3 ). The boundary conditions are r(0) = 0, 
r s (0) = xsin<9 + zcos<9, F(L) = 0, and M(L) =0 0. 
As in the lumped parameter model, v = uoz x r at steady 
state in the rod frame. 

The primary dimensionless groups governing the rod 
shape are the angle and the dimensionless rotation 
speed x = r)uL A /A = {LI if, where £ = [A/^uo)} 1 ^ 
is the characteristic length scale determined by bending 
resistance and viscous drag 0,0]. In addition, the aspect 
ratio L/a and the non-dimensional gravitational force 
gfi ro dL 3 /A are included in the analysis, but not para- 
metrically explored since they play a minor role. Fig- 
ure [3] shows that experimental measurements using two 
rod lengths collapse well onto a single curve when dimen- 
sionless speed, is plotted against dimensionless motor 
torque, M m L/A. The open symbols represent constant 
velocity rotation and trace out the entire S-shaped curve. 
The filament shape is stable at every prescribed value of 
X, and the shape changes continuously from the slightly 
bent shape to the helical shape as x increases. The closed 
symbols represent constant torque operation. There is a 
jump in rotation speed and filament shape at two dif- 
ferent torque values, depending if torque is ascending or 
descending. For descending torque, the time to reach 
steady state is prohibitively long, and the diamond sym- 
bols above the curve in Fig. [3] represent shapes that are 
still relaxing slowly to steady state. 

It is interesting to note that the nonlinear behav- 
ior of the speed-torque curve displayed in Fig [3] can 
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FIG. 4: (Color online) Steady-state shapes of a rotating rod 
from experimental measurements for x — 1-38 (O)? 4.25 (□), 
5.91 (0) (before transition) , and 164.63 (v) (after transition) , 
along with the shapes calculated from the nonlinear (— ) the- 
ory. The rod has L = 210 mm and 6 = 20°. 

be qualitatively explained using linear approximations 
valid for small rod deflections, since the shape of the 
rod depends nonlinearly on x even when the rod de- 
flection is small. For small the rod is aligned mainly 
along the z-axis and, disregarding gravity, the deflection 
r±(z) = (x(z),y(z)) obeys 

-^+zxr 1 = 0. (6) 

The solution to Eq. (|6|) is a generalization of Machin's 
solution to the in-plane bending problem Q, and of the 
same form as the solution for a flexible rod held parallel 
to but some distance from the axis of rotation [171 ]. 

To calculate the motor torque M m required to rotate 
the rod at speed cj, observe that the moment due to vis- 
cous drag must equal the elastic moment at the base of 
the rod: 

M m = -z J r x f ds = -Az • r s x r ss (0). (7) 

The second equality of (|7j) follows from ([4HS]) • The re- 
sults of the linear calculation for driving torque vs. speed 
are shown in Fig. [3] along with the experimental data for 
= 26°. For small x? M m increases linearly with x- We 
may use a simple scaling argument to capture the de- 
pendence of M m on x for x > 1. For large x? Eq. (|6|) 
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FIG. 5: (Color online) Theoretical calculation of the dimen- 
sionless propulsion force F p as a function of dimensionless M m 
for = 20°, using the nonlinear equations for the ideal case of 
zero gravity. The arrows denote the transition for ascending 
(|) and descending torque (j). 



on the motor, from the shape of the rod using F p = 
z • Jfds = — z • F(0). The kinematic reversibility of 
Stokes flow implies that a rigid rod rotating along the 
surface of a cone generates zero propulsive thrust. For 
small x, the elastic rod deforms slightly and generates lit- 
tle thrust. Above the critical torque, as the helical shape 
develops, the thrust increases abruptly (Fig. [5j). Since 
the shape of an actuated elastic filament cannot be de- 
coupled from swimming kinematics [20[, it would be an 
interesting generalization of our work to build an artifi- 
cial swimmer driven by a rotating elastic rod, tilted at 
the base to the rotation axis. 

This work was supported in part by National Science 
Foundation Grants Nos. CTS-0508349 (KSB) and NIRT- 
0404031 (TRP). As we were completing this letter, we 
learned of the work of the group of Fermigier, who have 
conducted a similar experiment [2lj |. 



implies that the shape of the rod is helical with an enve- 
lope that decays exponentially with length scale i. As- 
suming isotropy x ~ y and using force balance (|6j), the 
viscous force per length / ~ rjcoy ~ Ay/£ 4 , which im- 
plies a total viscous moment M V L/A ~ Ly 2 /£ 3 . On the 
other hand, the bending moment at the base of the rod 
scales as M^L/A ~ yL/£ 2 . Equating M v and Mb yields 
y ^ £ = Lx -1 / 4 0- Thus, for x ^> 1, the motor torque 
must scale as M m L/A ~ X 1 ^ 4 - Unfortunately, our exper- 
iment cannot access this high-speed regime due to limi- 
tations in the torque-speed characteristic of our motor. 
We also observe an intermediate scaling, M m L/A ~ x 1 ^ 2 - 
This scaling arises since in this sub-asymptotic regime the 
deflection y ~ Lx -1 ^ 4 , but the scale for bending of the 
rod is still L and not £ . 

For large the deflection of the rod is significant even 
for small x? an d the linear theory is inaccurate. How- 
ever, the nonlinear equations ([2H5]) are readily solved 
with shooting methods jTj|. The nonlinear theory gives 
a more accurate prediction for the speed-torque relation- 
ship in the high-speed regime where the linear and non- 
linear theories differ (Fig. [3]). As increases, the general 
appearance of the torque-speed relationship remains un- 
changed although both the critical torque and the jump 
in speed at the transition increase (not shown here). Fi- 
nally, for large x, the scaling analysis presented above re- 
mains valid, and the moment scales like x 1 ^ 4 for x ^> 1, 
just as in the linear theory. 

Figure [H shows the steady-state rod shapes for four 
different values of x = (£/^) 4 ? comparing experimental 
data (symbols) and the nonlinear theory (solid line). The 
agreement between theory and experiment is good. Note 
that the y-z projection shows how y(L) first increases 
with x and then decreases, in accord with our intuitive 
argument. 

We calculate the thrust, or axial force from the rod 
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